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Abstract. We have determined extinction corrections for a sample of 441 late-type stars in the inner Galaxy, which 
we previously searched for SiO maser emission, using the 2MASS near-infrared photometry of the surrounding 
stars. From this, the near-infrared extinction law is found to be approximated by a power law Ax oc y-i'9±o i 
Near- and mid-infrared colour-colour properties of known Mira stars are reviewed. From the distribution of the 
dereddened infrared colours of the SiO target stars we infer mass-loss rates between 10“^ and 10“® Mq yr”*^. 
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1. Introduction 

We study interstellar extinction in the direction a sam¬ 
ple of evolved late-type stars in the inner Galaxy (—4° < 
I < -1-30°, |6| < 1°), which we previously s earched for SiO 
mase r emission (“SiO targets” hereafter; iMessineo et al.l 
120021 Paper I). Since SiO maser emission reveals stellar 
line of sight velocities with an accuracy of a few kms“^, 
SiO maser stars are ideal for Galactic kinematics studies. 
Furthermore, the combination of stellar kinematics and 
physical properties of the stars, e.g. their intrinsic colours 
and bolometric magnitudes, will enable a revised kine¬ 
matic study of the inner Galaxy to reveal which Galactic 

Send offprint requests to\ M. Messineo, e-mail: 
mmessineSeso.org 

* This is paper no. 22 in a refereed journal based on data from 
the ISOGAL project. Based on observations with ISO (an ESA 
project with instruments funded by ESA member states and 
with the participations of ISAS and NASA), with telescopes 
of the European Southern Observatory in La Silla, Chile, and 
with the IRAM 30-m telescope on Pico Veleta, Spain. IRAM 
is supported by INSU/CNRS (France), MPG (Germany), and 
IGN (Spain). 

** Table 3 is only available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http: / / cdsweb.u-strasbg.fr / cgi-bin / qcat? J / A-\-A/ 


component and which epoch of Galactic star formation 
the SiO stars are tracing. 

A proper correction for interstellar extinction is of pri¬ 
mary importance for a photometric study of stellar popu¬ 
lations in the inner Galaxy, where extinction can be signif¬ 
icant even at infrared wavelengths. Interstellar extinction 
hampers an accurate determination of the stellar in trinsic 
colours and bolometric magnitudes I Messine^l2004ll . This 
is especially critical in the central Bulge region where in¬ 
terstellar extinction may exceed 30 visual magnitudes and 
due to the current uncertainty in the near-infrared extinc¬ 
tion law (30%) the uncertainty in the bolometric luminosi¬ 
ties of evolved late-type stars is at least 1 magnitude. 

The available near- and mid-infrar ed photometry 
of th e SiO t argets from the DENIS lEpchtein et al 


or th e biU t argets trom tne UFiNib llFp 
1994 |1. 2MASS llGutn et a,lJl200,‘^l . ISO GAlIc 


,_ _.Omont et al 

200,4 ISchuller et al.ll200.3lL and MSX l|Egan et al .I ll999t 
Price_et_^ .1 I 2 OOIII _™rveys were already presented by 


Messineo et al.l l|2004 Paper II). The sample consists 


mainly of large-amplitude variable (LAV) A symptotic 
Giant Branch (AGB) stars jPaperlll : IPanerJl . The esti¬ 
mates of interstellar extinction toward this class of objects 
are complicated by the presence of a circumstellar enve¬ 
lope with variable optical depth. Therefore, in order to dis¬ 
entangle circumstellar and interstellar extinction for each 
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AGB star one also needs to study the dust distribution 
along its line of sight. 

In Sect. [21 we discuss the uncertainty of the extinction 
law at near- and mid-infrared wavelengths and the con¬ 
sequent uncertainty of the stellar luminosities. In Sect. 0 
we describe the near-infrared colour-magnitude diagrams 
(CMDs) of field stars (mainly giants) toward the inner 
Galaxy and use the latter to derive the median extinction 
in the direction of each target. In Sect.0]we review the lo¬ 
cation of known Mira stars on the CMDs and colour-colour 
diagrams. In Sects. O and Owe use the median extinction 
from surrounding field stars to deredden our SiO targets 
and discuss their colours and mass-loss rates. The main 
conclusions are given in Sect.O 

2. Interstellar extinction law 

The composition and abundance of interstellar dust and 
its detailed extinction properties remain unclear, thus lim¬ 
iting the accuracy of stellar population studies in the inner 
Galaxy. In the following we discuss the near- and mid- 
infrared extinction law, in order to assess the uncertainty 
in the extinction correction. 



5 10 15 20 25 
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Fig. 1. Filter transmission curves and extinction laws 
as function of wavelength. The continuous line (Gurve 1 ) 
shows a fit to the values (dots) given bv iMathi^ lll99(1l : 
th e dashed curve shows t he parametric expression given 
bv [R osenthal et"^ 1 li2000i) plotted using a value of the sil¬ 
icate peak 219 , 7 /^ 2.12 = 1.0 (Gurve 2). The latter is also 
shown wi thout the m inimum around 4-8 ^m (Curve 3), 
following iLut j (^9^. In the top panel t he transmis¬ 
sion c urves of the MSX A,C,D and E filters llPrice et alJ 
l2nniii are shown, an d in the bottom pa nel those of 
the ISOCAM filters llR1omma,ert et a,f]l2flfl3ll used in the 
ISO GAL survey. 


2.1. Near-infrared interstellar extinction 


Interstellar extinction at near-infrared wavelengths (1- 
5 ^m) is dominated by graphite grains. Near-infrared 
photometric studies have shown that the wavelength- 
dependence of the extinction may be expressed by a 
power law A A (x where o wa s found t o range be¬ 
tween 1.6 (|Rieke_fcLebofsk3[l9 ^ an d 1.9 l|Glasslll9^ 
iLandini et al.lll984Hvan de Hulstlll94(jl . 

One can estimate the near-infrared extinction by mea¬ 
suring the near-infrared reddening of stars of known 
colour. When using broad-band photometric measure¬ 
ments, one needs to properly account for the bandpass, 
stellar spectral shape, and the wavelength-dependence of 
the extinction. As a function of the band extinction we 
therefore computed a grid of “effective extinction” values 
for the DENIS I and 2MASS J, H, and bandpass. 
This effective extinction was computed by reddening an 
MO III stellar spectrum llFluks et al.lll 994^ with a power 
law extinction curve and integrating over the respective 
filter transmission curves. When we convolve the filter re¬ 
sponse with a stellar sub-type spectrum different from the 
MO III, only the effective I-band extinction values change 
significa ntly, e.g. decreasing by 3% for a M7 III spectrum 
(see also Ivan Loon et ^l2003|) . 

The Ag-band extinction Ak^ can then be expressed as 


Aks = CjK X E{J — Ks), 


Aks = Chk X E{H - Ks), 

where E{J — Ks) and E{H — Ks) are the reddening in 
the J — Ks and H — Kg colour, respectively, and the 
C are constants. These relations are independent of vi¬ 
sual extinction and of the coefficient of selective extinc¬ 
tion, i?v = Ay/E{B — V), but they depend on the slope 
of the near-infrared extinction power law (see Table ^ . 
However, to provide the reader with the familiar ratio 
between near-infrared effective extinction and visual ex- 
tinctio n, we also used the c ommonly adopted extinction 
law of iGardelli' et al.l lll989ll . Such ratios may be useful 
in low-extinction Bulge windows, where visual data are 
also available. The analytic expression given by Cardelli 
et ah, which depends only on the parameter Ry, is based 
on multi-wavelength stellar colour excess measurements 
from the violet to 0.9/rm, and extrapolates to the near - 
infrared using the power law of iR.ieke fc Lebofskvl (1198511 . 
We extrapolated Gardelli’s extinction law to near-infrared 
wavelengths using a set of different power laws (Table P) . 

Uncertainty in the slope of the extinction law produces 
an uncertainty in the estimates of the near-infrared extinc¬ 
tion of typically 30% in magnitude (see Table P). For a Ks 
band extinction of Ak 3 = 3 mag the uncertainty may be 
up to 0.9 mag, which translates into an uncertainty in the 
stellar bolometric magnitudes of the same magnitude. 

In Sect. 121 we show that a power law index a = 1.6 
is inconsistent with the observed near-IR colours of field 
giant stars toward the inner Galaxy, and that the most 
likely value of a is 1.9 ± 0.1. 
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Table 1. Near-infrared effective extinction, (A)/Ay oc A for various filters. A different value of Ry = Ay/E{B-V) 
does not affect Aks/E(H — K^) and Ak,,/E(H — K^), but the slope of the extinction curve does. Our findings favour a 
model with a = 1.9 (see Sect. 0). 


Ai/Av 

Aj/Av 

Ah/Av 

Aks / Av 

Ak3/E(J-KJ 

AkJE(H-KJ 

a 

Rv 

Ref. 

0.592 

0.256 

0.150 

0.089 

0.533 

1.459 

1.85 


Glass 119991 

0.482 

0.282 

0.175 

0.112 

0.659 

1.778 

1.61 

3.09 

Rieke Lebofskv l'l^<j)(>l 

0.606 

0.287 

0.182 

0.118 

0.696 

1.842 

1.61 

3.10 

Card el li et al. ("^22^ 

0.563 

0.259 

0.164 

0.106 

0.696 

1.842 

1.61 

2.50* 

5 ? 

0.606 

0.267 

0.158 

0.096 

0.561 

1.548 

1.85 

3.10 

Cardelli et al. 119891 + 

0.563 

0.240 

0.142 

0.086 

0.561 

1.548 

1.85 

2.50 

?) 

0.606 

0.263 

0.153 

0.092 

0.537 

1.496 

1.90 

3.10 

5 ) 

0.563 

0.237 

0.138 

0.083 

0.537 

1.496 

1.90 

2.50 

5 ) 

0.606 

0.255 

0.144 

0.084 

0.493 

1.401 

2.00 

3.10 

5 ) 

0.563 

0.229 

0.130 

0.076 

0.494 

1.400 

2.00 

2.50 

?) 

0.606 

0.238 

0.127 

0.070 

0.420 

1.236 

2.20 

3.10 

?) 

0.563 

0.213 

0.114 

0.063 

0.420 

1.236 

2.20 

2.50 

5 ) 


* recent determination toward the Bnlge (e.g. IUdalskill20n3ll . 

■*" parametric expression modified to extrapolate to A > 0.9 fim with A““ 


2.2. Mid-infrared interstellar extinction 

The mid-infrared extinction (5-25 /im) is characterised by 
the 9.7 and 18 /rm silicate features. The strength and pro¬ 
file of these features are uncertain and appear to vary from 
one line of sight to another. A standard graphite-silicate 
mix predicts a minimum in A;^/A 2 .i 2 at 7/im, which has 
however been not ob served to be very pronounced toward 
the Galactic Centre l|lait7lll9^llnitz et a.iJIl99(i|l . 

The mid-infrared extinction curve of lMathi'3 ()l 99fl|l is 
commonly used with its combinatio n of a power l a w and 
the astronomical silicate profile of IPraine fc Tei (11984 
with A 9 . 7 /A 2.2 — 0.54 - a value found in the di ffuse in¬ 
terste llar medium toward Wolf-Rayet stars (e.g. iMathij 
Howe ver, from ob servations of hydrogen recombi¬ 
nation lines, llaitzl ll1999^ concluded that A 9 . 7 /A 2.2 — 1-0 
in the direction of the Galactic centre; and analysing 
the observed Ho level populations toward Orion OMC-1, 
iRosenthal et all ll2004 derived A 9 . 7 /A 2.12 = 1-35. 

By using a parametri c mid-infrared extinction curve 
bv lRosentha.1 et all 1 200fl(l . we computed two curves with 
A 9 . 7 /A 2.2 = 1 . 0 , one in combination with the minim um at 
4-8 um (Curve 2) and one without it as suggested bv lLiitd 
1 I 1999 I 1 (Curve 3) (see Fig.0. Using these extinction laws 
we computed the effecti ve extinctions (A) / Aks for the M- 
type synthetic spectra IIFluks et a1.l 11994 in the various 
filters used by the ISOGAL and MSX surveys (Table 0. 
The ratios (A) /Ak^ are not sensitive to the stellar sub- 
type used. An increase in the ratio 49 , 7 / 42.2 from 0.54 to 
1.0 results in an increase between 0.15 and 0.20 x4k 3 of 
the average attenuation in the LW3, LW9, C, D, and E 
spectral bands. The spectral bands of the LW2 and LW5 
filters are not very sensitive to the intensity of the silicate 
feature, but to the minimum of the extinction curve in the 
4-8 fim region. Although {A)/Ak^ varies with Ak^, these 
variations are small compared with those arising from dif¬ 
ferent choices of the mid-infrared extinction law. 


In the most obscured regions (4 k„= 3) uncertainties 
for the ISOGAL and MSX filters range from 0.45 mag 
{LW3, LW9, D) to 0.85 mag (4). Because their energy 
is emitted mostly at near-infrared wavelengths, this has 
a negligible effect (0.1 mag in average) on the calculated 
Mboi of the SiO targets and, therefore, also on the mass- 
loss rate estimates (see Sect. 0. In the following we use the 
Lutz law (Curve 3), since it ensures consiste ncy between 
mid-infrared and near-infrared stellar colours llJiang et alJ 

3 . Interstellar extinction of field stars derived 
from near-infrared colour-magnitude diagrams 

Most of the sources detected by DENIS and 2MASS to¬ 
ward the inner Galaxy are red giants and AGB stars. 
Because the intrinsic colours of giants are well known and 
steadily become redder with increasing luminosity, one can 
study the near-infrared GMDs to estimate the average ex¬ 
tinction toward a given line of sight for a population of 
such stars. This approach has been used to map the ex- 
tinction^in^the^cmUal^region^of^the^GM axy {\l\ < 10°) 

Assuming the SiO target stars to be spatially mixed 
with the red giant stars, we can similarly estimate extinc¬ 
tion, Aks, toward the 441 SiO targets. In both the (Ag, 
J—As) and (Ag, A—Ag) diagrams the colour and magni¬ 
tude of field stars within a 2-4' radius from each SiO target 
were shifted to the reference red giant branch, RGB, (see 
Sect, mil along the reddening vector. Selection effects due 
to magnitude limits were taken into account ('Sect. I!01l . 
Median extinction of the field was thereby determined in 
both the (Ag, J—Ag) and (Ag, A—Ag) planes. To exclude 
foreground stars, an iterati ve 2a clipping wa s applied to 
the extinction distribution dPutra, et al.ll2flfl^ . 

Nearly 90% of the SiO target stars appear redder than 
their neighbouring field stars (Fig.0, which implies that 
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Table 2. Effective extinction, {A)/Aks, using M-giant spectra llFluks et al iim for different bands defined by the 
ISOCAM and MSX filters (see Fig.^J. AksI-^ 2.12 = 0.97. 


Filter 

Aref 

fim 

AA 

fim 

Curve 1 (Mathis) 
{Ag. 7 /A2.12 ~ 0.54) 
{A)/Aks 

Curve 2 

{Aq. 7 /A2.12 = 1.00) 
{A)/Aks 

Curve 3 (Lutz) 

(Ag, 7 /A 2 .i 2 = 1.00 & no minimum) 
(A)/Aks 

LW2 

6.7 

3.5 

0.21 

0.21 

0.41 

LW5 

6.8 

0.5 

0.18 

0.15 

0.41 

LW6 

7.7 

1.5 

0.21 

0.26 

0.43 

LW3 

14.3 

6.0 

0.18 

0.34 

0.34 

LW9 

14.9 

2.0 

0.14 

0.29 

0.29 

A 

8.28 

4.0 

0.26 

0.38 

0.55 

C 

12.1 

2.1 

0.25 

0.49 

0.49 

D 

14.6 

2.4 

0.14 

0.29 

0.29 

E 

21.3 

6.9 

0.17 

0.41 

0.41 


the SiO stars are intrinsically obscured, provided their 
spatial distribution is the same as that of RGB stars. 

3.1. Reference red giant branch 

Our method for measuring the extinction needs reference 
to an RGB: we adopted that of 47 Tuc, bro ught to 8 kpc 
by ad opting a distance modulus of 13.32 llFerraro et alJ 
[1991). This choice is supported by a recent new near- 
infrared col our-magnitud e analy sis of the RGB in globular 
clusters bv I Valenti et alJ 1 I 2 OO 4 I 1 . who obtained {J,H,K) 
photometry taken at the ESO 2.2 m telescope with a spa¬ 
tial resolution of 0.3 and 0.5" pixel“^, half that of 2MASS. 
They showed that within 0.05 mag the RGB of 47 Tuc and 
the metal-rich Bulge globular clusters (such as NGC 6528, 
NGC 6553, NGC 6540) have identical colours. 

In the 2MASS Kgo, (J — Ks)o CMD, the upper part 
{Ks< 12 mag) of the 47 Tuc RGB is well represented by 
a linear fit 

(J - K,)o = 2.19(±0.02) - 0.125(±0.002)Kso; (1) 

in the {H — Ks)o vs. Kso diagram a second-order polyno¬ 
mial fits well the RGB {Kso < 12 mag): 

(H - K^o = 1.73(±0.22) - 0.268(±0.05) x Kso (2) 

-f0.011(±0.002) X kIq. 

The Ks, {H — Ks) plane is less sensitive to extinction 
than the Kg, (J — Kg) plane: while a colour change of 0.1 
in (J — Ks)q implies a change of 0.05 in Ak^ (0.6 in Ay), 
a change of 0.1 mag in {H — Kg)Q implies a change in 
0.15. A shift in distance modulus of the reference 
RGB of ±2 mag results in a small change in the extinction 
of AksTO-IS mag. 

3.2. Determination of extinction value and extinction 
law in the J, H, Ks CMD 

To properly estimate the interstellar extinction along a 
given line of sight one should consider only the region of 
the CMD where the upper RGB is well defined, i.e., not 
affected by large photometric errors or by the diagonal 


Table 3. Extinction values. The identificati on numb er 
(ID) of the SiO target, as in Table 2 and 3 of lPanerl is 
followed by the field extinction Ak^ , by the corresponding 
dispersion of individual extinctions of field stars and by 
the total (circumstellar plus interstellar) extinction in Kg- 
band toward the target star (tot). The flag (Eg) is unity 
when the SiO target is classified as a “foreground object”. 


ID 

mag 

'^AKs 

mag 

tot 

mag 

Fg 

1 

0.96 

0.18 

1.63 


2 

1.23 

0.31 

2.81 


3 

1.90 

0.33 

1.79 


4 

2.16 

0.52 

3.00 


5 

1.51 

0.25 

1.75 


6 

0.59 

0.09 

1.16 


7 

0.57 

0.07 

0.51 


8 

1.98 

0.41 

1.53 

1 

9 

2.14 

0.32 

2.59 


10 

0.61 

0.11 

1.48 



* The full table is available in electronic form at 
the CDS via anonymous ftp to cdsarc.u-strasbg.fr 
(130.79.128.5) or via http : //cdsweb.u — strasbg.fr/cgi — 
bin/qcat?J/A -|- A/(vol)/(page). 


cut-off from the 2MASS detection limits (see Fig. [2Jl that 
would bias the calculatio n of the median extinction toward 
a lower value (see als o ICotera et all l2n0(Tt IPutra et alJ 
120031: iFiger et al ] |2004ll . 

In the inner Galaxy, the photometric error is typi¬ 
cally smaller than 0.04 for stars with Kg < 12 mag. To 
quantify the incompleteness due to the diagonal cut-off 
with zero extinction, our average 2MASS detection lim¬ 
its of J = 16.0 and H = 14.0 correspond to a RGB Kg 
magnitude of 15.2 and 13.0 mag, respectively, at a dis¬ 
tance of 8 kpc. Accounting for scatter in the observed 
colours of ±0.5 mag, the RGB would be well sampled to 
J = 15.5 and H = 13.5 mag, corresponding to Ks< 14.6 
and 12.6, respectively. With a Kg extinction of 3 mag 
(~ 5.0 mag in H, ~ 8.6 mag in J), which is a typ¬ 
ical value toward the Galactic centre, these RGB com- 
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#178, #179. #181, <I =9.8‘, b = 0.1*) 


1.475 
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4 6 a 100 

J-K* 
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Fig. 2. Colour-magnitude diagrams of 2MASS datapoints (small dots) of good quality (0 < flag-red < 3) located 
within 4' from the position of the SiO target (big dot). Three fields at equal median extinction are combined in each 
panel. The right-hand continuous line indicates the locus of the reference RGB curve (see Sect. EU, adopting a distance 
of 8 kpc and reddening it with the median extinction of field stars (Ak^)- The left-hand continuous curve shows the 
tracels of clump stars for increasing distance an d extinction along a given line of sig ht (see Sect. JO, obtai ned using 
the extinction model bv IPrimmel et a.H 1 200,311 and the absolute magnitudes from IWa.inscoat et al.l ill 992ll . Dashed 
lines indicate the diagonal cut-off due to the detection limits in J and H. Circles on the clump trace mark a distance 
from 1 to 5 kpc with a step of 1 kpc downward. 


pleteness limits would rise to Ks = 7.1 and 10.8 mag in 
the (itTg, J — Ks) and {Kg, H — Ks) planes, respectively. 
Therefore only the H band can provide a sufficient num¬ 
ber of red giant stars to match the reference RGB. With 
the 2MASS data {Kg, J — Kg) CMDs are useful for extinc¬ 


tion determinations only to a Kg extinction of about 1.6 
mag. 

For fields with Ak 3 < 1-6 mag Table El lists the ex¬ 
tinction values determined from the {Kg, J—Kg) plane, at 
larger extinction values those from the {Kg,H—Kg) plane. 
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0 0.5 1 1.5 2 8.5 


Ak (H.KJ 

Fig. 3. Comparison between the median extinction ob¬ 
tained from the {Ks,J—Ks) and {Ks,H—Ks) planes. Only 
sources detected in J,H,Ks above the Kg completeness 
limits are used. In each panel a different value of the spec¬ 
tral index of the extinction power law, a, is adopted. 


By using both the {Ks,J—Kg) and {Ks,H—Ks) planes 
one can in theory fit the observed field star colours to the 
RGB and ht both the absolute average Kg band extinction 
and the index, a, of the extinction power law. For this 
we selected field stars which were detected in J, H, and 
TOand which are brighter than the Kg completeness limits 
for the RGB at the extinction of each held. For each held 
we determined the median of the extinctions of the held 
stars in both the {Ks,J — Kg) and in the {Kg,H — Kg) 
CMDs. The power law index a was chosen such that the 
extinctions agreed in the two planes. For diherent values 
of a. Fig. 13 shows the differences in median extinction 
values AK^{J,Kg)—AK^{H,Kg) plotted against Ak^. Only 
for a = 1.9 ± 0.1 are the two extinction estimates found 
to be consistent. 

In the followi ng we use a = 1. 9, which is co n sisten t 
with the work of IClas i (l^ 999^ and ^Mdm^^a I] lll984ll . 
and with the historical Curve 15 of van de Hulsfj lll946h . 

3.3. Outside the Bulge 

K2 giant stars are the dominant populati on of late-type 
stars seen along the galactic disk (e.g. IPrimmel et al.1 
l2003t iLopez-Corredoira et alJlio?)^ . They correspond to 
red clump stars in metal-rich globular clusters such as 47 
Tuc. The location of clump stars on the CMD is marked 
in Fig. [3 


Toward the Bulge, the Bulge RGB population is dom¬ 
inant; therefore, the median interstellar extinction is not 
affected by possible foreground clump stars. This is not 
the case in the disk {I 10°), where one must ignore the 
foreground clump stars before fitting the RGB. As likely 
clump stars we identify those within 0.3 mag from the 
J—Kg colour of the clump trace, and as giants those stars 
redde r than the clump stars (e.g. iLonez-Corredoira et alJ 

I 2 OO 2 I) . 

3.4. Dispersion of the extinctions along a line of sight 

Along with the median extinction of the field stars, Ak^, 
we also determine the standard deviation of the distri¬ 
bution of the individual extinction estimates, ■ The 
patchy nature of the extinction is apparent even within the 
2—4' radius sampling area. This patchiness integrated over 
a longer path generates larger ctak^, with increasing extinc¬ 
tion for Bulge lines of sight. The Icaks uncertainty in the 
field extinction varies from ~ 0.2 mag when Aks= 0.6 
mag up to ^ 0.7 mag in the regions with the largest ex¬ 
tinction (Ak„> 2.0). In fields at longitudes longer than 
10°, a larger ctaks i® found than in Bulge helds of similar 
median extinction. This is probably due to the presence 
of several Galactic components, e.g. the disk, arms, bar, 
and molecular ring. 


4. Near-infrared properties of known Mira stars 


Presently the pulsation periods and amplitudes of our 
SiO target s are un known, although most of them must 
be LAVs ( Paper Illl . This is suggested by photometric 
variations seen in measurements taken at different times 
(DENIS. 2MASS and MSX) by the strong 15 »m em ission 
( Paper ijl and by the SiO maser emission ij PaperJl . 

Although they are 20 ti mes less numerou s than 
semiregular AGB stars (SRs) ijAlard et af.ll2nni|i . Mira 
stars are among the best studied pulsating variable stars. 
They are regular long period variables (LPV) with visual 
amplitudes larger than 2.5 mag and K band amplitudes 
exceeding ~ 0.3 mag (LAV stars). Since large amplitudes 
tend to be associa ted with the most regular light curves 
l|Cioni et af .11200,'^ . amplitude is the best choice for clas¬ 
sifying Mira stars. 

To analyse the colours of our SiO targets, in partic¬ 
ular to check the quality of the extinction corrections, 
it is useful to have comparison samples of well studied 
large amplitude LPV stars with a wide range of intrin¬ 
sic colours, at low interstellar extinction and at differ¬ 
ent Galactic positions. Therefore, we e xamined: two sam¬ 
ples of Mira stars in th e solar vicinity ifOlivier et a,1.ll2fl(Til 
IWhitelock et al.l l2000l): one samp l e towa rd the Galactic 
Gap taken from IWhitelock et a.1 1 (11994 ; 18 Mira stars 
dctactaci by T ]i 99^ i n the Sgr-I field; 

104 IRAS Mira stars in the outer Bulge dWhitelock et alJ 
Em 3); and a sample of LPV in the Large Magellanic 
Gloud l|Whitelock et ahlliool) . All these stars have IRAS 
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12 /rm magnitude [12] and mea n J, H, K magnitudes in 
the SAAO system ilCarteilll99(il . 

Stellar fluxes a re given already corr ected for reddening 
only in the w ork oflO livier et alJ 1120 0111. For the L PV stars 
analysed by IWhitelock et all l|2000l 1 199"^ l2003ll the ef¬ 
fects of interstellar extinction are negligible because these 
stars are nearby or outside of the Galactic plane, so we 
did not correct these for extinction. We de reddened the 
Baade Sgr-I window data llGlass et al ■ll 99,^ adopting our 
favourite extinction curve {a = 1.9) a nd Ak. = 0.15 mag , 
consistently with the value adopted in iGlass et alJ l|l991ill . 
We corrected the mag nitudes of the outer Bulge Miras 
llWhitelock et al.lll99lli for reddening adopting values of 
Aks derived from their surrounding stars (see next Sect.). 
Next we analyse the location of these well-known Mira 
stars in the near-infrared CMDs and colour-colour dia¬ 
grams. 

4.1. Colour-magnitude diagram of outer Bulge Mira 
stars and surrounding field stars 



02468 01234 




Fig. 4. Dereddened colour-magnitude diagrams. Big 
dots represent the o uter Bulge Mira stars found by 
IWhitelock et alJ lll99 jl . The magnitudes plotted are mean 
magnitudes at the equal distance o f 8 kpc, adopting the 
distances of lWhitelock et ahl 1 199l|l . Small dots represent 
the point sources detected by 2MASS within 1' from each 
Mira star. 


Th e outer Bulge LPV stars studied bv lWhitelock et aP 
are particularly interesting for a comparison with 
our SiO targets because they were selected on the ba¬ 
sis of their IRAS fluxes and colours with crit eria simi¬ 
lar to t hose used to select our MSX targets llPaoer ij : 
IPaDerJ] . Since their main periods range from 170 to 
722 days and their K amplitudes from 0.4 to 2.7 mag, 
they are classical Mira stars. Their distribution of dis¬ 
tance moduli, estimated from the period-luminosity rela¬ 


tion ilWhitelock et al.lll99lll . peaks at 14.7 mag with a 
a 0.5 mag. Since they are at latitudes between 6 and 
7°, they are in regions of low interstellar extinction. All 
this makes them ideal for a comparison with our SiO tar¬ 
gets, the study of which is complicated by the large inter¬ 
stellar extinction at their low latitudes. 

An RGB is clearly apparent in Fig. 0] which shows 
the 2MASS point sources within 1' of each Mira star. The 
median extinction toward each held ranges from AKg = 
0.01 to 0.30 mag with a typical dispersion of 0.01-0.08 
mag. On the CMDs the Mira stars appear mostly brighter 
than the RGB tip of t he held stars (K = 8.2 mag at a 
distance of 8 kpc, see iFrogel fc WhitfordlllQj^ . Due to 
the presence of a circumstellar envelope, Mira stars have 
red colours (up to {H — Ks)o = 3 mag) and lie on the red- 
side of the giant sequence. It is therefore not possible to 
derive the interstellar extinction toward Mira stars from 
their colours relative to the RGB. 

There is no reason to assume the Mira stars are spa¬ 
tially distributed differently from other giant stars, so 
that the extinction toward the surrounding field stars may 
serve as an approximation of that of the respective Mira 
star. One worry in this assumption is not knowing the 
actual distribution of extinction along the line of sight. 



Fig- 5. D ereddened colours of known Mira stars 
jGlass^^^ ] ll99^ lOlivier et al.1 l200lt IWhitelock et alJ 
Ii99iL l2nnni l^994^ . Near-infrared mean magnitudes are 
used. The two curves represent M3 (upper) and MIO 
(lower) type stars with increasing mass-loss rates (indi¬ 
cated by squares and labels xl0 ~ ^ M^Tj /vear). as mod¬ 
elled by iGroenewegen fc de ,lond ill 99311 . The cross in¬ 
dicates the position of an MIO star without mass-loss 
{{H — K)o = 0.48 and {J — K)q = 1.34). The arrow shows 
the reddening vector for Ak= 1 mag. 
























































M. Messineo et al.: 86 GHz SiO maser survey of late-type stars in the Inner Galaxy 


4.2. Colour-colour diagram of Mira stars 



Fig. 6. Dereddened colours of infrared-monitored Mira 
stars (based on the IRAS 12 gm and near-infrared 
mean magnitudes): in the solar vicinity (filled squares) 


JOlivier et al. booill: detected by Hie 

parcos (open trian- 
the South Galactic 

fflesl JWhitelock et al.l 200fil: toward 

Gap (open squares') JWhitelock et al.1 

199j); in the Baade 

Sgr-I window (biff dots') JGlass et al. 

in the outer 

Bulffe Ismail dots') d’Whitelock et all 

199lll; in the Larffe 

Maffellanic Gloud (’crosses') l|’\Yhif^l 9 (;k 9,lJll994|). The 


dotted line is th e best fit to an IRAS sample of oxygen- 
rich AGB stars llvan Loon et al ] |l997ll . The arrow shows 
the reddening vector for Ak= 1 mag. 


For Mira stars with low mass-loss rate (< 10 ' Mq 
vr~ ^). the (J — K)r[ co lours range between 1.2 and 


1.6 dWhitelock et al.l[200flll . Dust-enshrouded IRAS AGB 
stars with mass-lo ss rates of 10“® — 10“^ Mq yr“^ 
are much redder, with { J — K)q rang- 


ijOlivier et al.ll2f 
ing from 2 to 6.5 mag. The overall distribution (Fig. EJ 
appears to form a sequence of ever redder colours with 
increasing mass-loss rate, a trend that is well reproduced, 
e.g., by a model for an MIO type AGB star with increas¬ 
ing shell opacity I Groenewegen fc de Jondll993ll . Thus, a 
higher mass loss has the same effect on {J — K)q and {H — 
K)q colours as more interstellar absorption/reddening, 
thus making a distinction between intrinsic and interstel¬ 
lar reddening in these colours impossible. 

In contrast, in the (J — K)q vs. {K — [I2])o and the 
{H — K)o vs. (AT — [12])o diagrams (Figs. El andCj), a sep¬ 
aration between interstellar and circumstellar extinction 
could theoretically be made. However, this is hampered by 
the large dispersion of the Mira stars around the colour- 
colour fiducial sequence (0.5 mag), which is due to the 
non-contemporaneity of the near- and mid-infrared data 
and to the dependence of such relations on metallicity 
and stellar spectral type. The uncertainty of interstellar 



(H-K)„ 


Fig. 7. Dereddened colours of infrared-monitored Mira 
stars, based on the IRAS 12 /im and near-infrared mean 
magnitudes. Symbols are the same as in Fig. El The dot¬ 
ted line is t he best fit to an IRA S sample of oxygen-rich 
AGB stars llvan Loon et al ] ll998ll . The continuous curve 
is our best fit for Galactic Mira stars. The arrow shows 
the reddening vector for Ak= 1 mag. 


extinction estimates found by shifting a photometry point 
along the reddening vector onto the H — K,K — [12] curve 
is larger than Ak= 1 mag for K — [12] > 3.5 mag. 

In the {H — K)q vs. {K— [12])o plane known Mira stars 
appear to lie below the fiducial colou r-colour sequence of 
IRAS-selected oxygen- rich AGB st ars llFouaue et al.ll992t 
iGuglieliiio et alJll9^ derived bv Ivan LoonendTllbb^ 
The off set could be due to the H band water abso rption 
bands I Frogel fc Whitfor J Il987t iGlass et "^ Il995ll that 
are found to be strong in LAV stars. However, we did not 
find a correlation between colour and variability index for 
these IRAS-selected stars. The offset might also suggest 
a problem transforming the original ESO photometry of 
such cold stars in the SAAO system. 

Mira stars with 0.2 < {H — K)q < 3 mag fit 


{K - [12])o = 4.26(±0.04) -k 5.95(±0.I7) log(i7 - K)o, 

with an rms d eviation of 0 .5 mag. Although we use the 
SAAO system llGarteill990ll . in Appendix A we show that 
this fit is also valid in the 2MASS photometric system. 


5. Interstellar extinction and intrinsic colours of 
the SiO targets 

Due to interstellar extinction along the line of sight, for a 
given Ag —[15], the H—Kg colours of the SiO targets are 
much redder than those expected from the colour-colour 
relation of known Mira stars (upper panels of Fig.E|). 

The distribution of the SiO targets in the dereddened 
{H — Ks)o, {K— [15])o diagram (lower panels of Fig.El) us¬ 
ing the “field” extinction values, approaches that of known 
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Fig. 8. Left upper panel: 2MASS J — H versus H—Kg colours. Stars with upper magnitude limits are not shown. 
Dots and triangles represent objects with Kg smaller and larger than 6.0 mag, respectively. Filled and open symbols 
represent SiO detections and non-detections, respectively. “Foreground objects” were removed. The arrow shows the 
reddening vector for Ak„ = 1 mag. Right upper panel: 2MASS Ag —[15] versus H—Kg colours. The curve represents 
the best fit to the colours of known Mira stars (see Sect. ^J.Left lower panel: Dereddened 2MASS (J — H)q versus 
{H—Ks)q colours. Right lower panel: Dereddened 2MASS (Ag—[15])o versus (A—Agjo colours. 


Mira stars. Therefore, the “field” median extinction is a 
good approximation of the interstellar extinction for most 
of the SiO targets. Though the dispersion of individual 
field star extinctions along a given line of sight is consid¬ 
erable (from 0.1 to O.SctakJj the distribution is strongly 
peaked, especially in the Bulge region. Furthermore, the 
lifetime of a star on the AGB evolutionary phase is very 
short: it is about 5% of the time spent on the helium core 
burning phase and from 0.1 to 2% of t he time spent on the 
maiii sequence phase of such a star llVassiliadis fc Wo^ 
llQO.ll) . We therefore expect most AGB stars to be located 
in the region with the highest stellar density along the line 
of sight. 


There is still an asymmetry in the distribution of the 
SiO targets around the fiducial colour line of known Miras, 
which suggests that we could have underestimated the in¬ 
terstellar extinction for part of the sample up to 15%. In 
regions of high extinction, Mira-like stars are detectable 
at larger distances than ordinary field stars due to their 
high near-infrared luminosity. Deeper infrared observa¬ 
tions ar e needed to obtai n more accurate extinction es¬ 
timates i Figer et al ] |2004 . 


5.1. “Foreground objects" 

We can measure the total (circumstellar plus interstellar) 
extinction of each SiO target (Table EJ, by assuming a 
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stellar photospheric (J — Ks)o colour of 1.4 mag and an 
{H — Ks)o colour of 0.5 mag (cf. Fig. O. 

The difference between these total extinction estimates 
and 2MASS and DENIS J, Kg data gives an rms of AAk^ 
= 0.2 mag, while estimates from 2MASS H—Ks and J—Kg 
colours give an rms difference of AAk = 0.13 mag. For 
only 12 SiO targets we do not have any observed H — Kg 
or J — Kg values to determine the total extinction. 

As expected, on average the target stars show larger 
total extinctions than their surrounding field stars. This 
is not the case, however, for a group of ^ 50 mostly 
very bright target stars {Kg < 6.0) at various longitudes, 
marked with flag Fg= 1 in Table El which have total ex¬ 
tinctions lower (a least Ictaks) than the “held” extinc¬ 
tion and are therefore likely to be foreground objects. We 
dereddened the “foreground objects” by directly shifting 
them on the iJ, Kg, [15] colour-colour sequence. 


6 . Intrinsic colours and mass-loss rates 


The stellar mass-loss rate is best estimated from mea¬ 
surements of CO rotational lines. The CO emission arises 
in the circumstellar shell. However, because of confu¬ 
sion with interstellar CO emission, it is difficult to ob- 
tain such measurement s toward stars in the inner Galaxy 
i Winnberg et al Although infrared emission also 

arises in the stellar photosphere, stellar outhows may be 
studied using the infrared emission of dust grains that 
form in the cool circumstellar envelopes. Relations be¬ 
tween the infrared colours (e.g. J — K, K — L, K — [12] 
or AT — [15]) of 0-rich AGB stars and t heir mass-loss rate 
have been established empirically (e.g. lAlard et al.ll200lt 
lOlivier et alJl200it IWhitelo(^ et al.lfl994l) and supported 

bA^theoretical^imdds (e.gJ^ro^^w^^nJ^^lg_Jo^L993t 

The empirical relation between the {K — [15])o colour 
and mass loss rate, M, is very useful to study stars 
detected in the 2MASS or DENIS surveys and in the 
ISOGAL or MSX surveys toward the most obscured re¬ 
gions of the Galaxy. 

The uncertainties arising from the variability of the 
stars and the temporal difference between the Kg and 
15/rm measurements are somewhat alleviated by using an 
average of the 2MASS and DENIS Kg fluxes and of the 
ISOGAL and MSX 15/J,m measurements. The remaining 
r.m.s. uncertainty of the mass-loss rate is t hus a factor ^2 
for M > 10"® Mq yr"! JOiha et al.ll2003ll. _ 

Following the prescription of l.leong et al.1 11200311 and 
lOiha et al.l 11200311 we obtained mass-loss rates for the SiO 
targets, the distribution of which is shown in Fig. M 90% of 
the sources have implied mass-loss rates between 10"^ and 
2 X 10"® Mq yr"^, with a peak in the range 10"® - 10"® 
Mq yr"^, although the apparent distribution is widened 
by the uncertainty of the mass-loss-colour relation, by the 
photometric uncertainty, and by the effects of variability. 

Our selection criterion on (Ag—[15])o for the ISOGAL 
sample completely eliminated sources with M > 10"®, 
while the elimination of OH/IR sources and the criteria 


on A — D and C — E colours for the MSX sample also 
considerably reduce d the proportion o f sources with large 
mass-loss rates (see IPauer lit IPaoer J) . The same results 
are obtained considering the subsample of targets with 
the best extinction corrections, i.e. with ctaks < 0-2 mag. 
Furthermore, the distribution of the mass-loss rates of the 
targets with detected SiO mase r emissio n appears similar 
to those without detected SiO llPaoer Jl . 

A possible underestimation of the interstellar extinc¬ 
tion results in an overestimation of the mass-loss rate. 
Adopting Mathis’ mid-infrared extinction law rather than 
Luts’s, the distribution of the mass-loss rates only slightly 
shifts toward lower values (see Fig. El). 



log(Mass loss rate) [Mq yr 


Fig. 9. Distribution of mass-loss rates derived from the 
A'—[15] vs. M relation fleong et ^ I 2003t . The continuous 
line shows the distribution for all SiO targets dereddened 
using Lutz’s extinction law (Curve 3) and the dotted line 
that of foreground stars. The dashed line is the distribu¬ 
tion of all SiO targets using the mid-infrared extinction 
law of Mathis (Curve 1). 


7. Conclusion 

We estimated the interstellar extinction toward 441 SiO 
maser AGB stars. From the 2MASS stars within 2-4' of 
each SiO star, we estimated a mean extinction toward 
the SiO stars. We therefore shifted the colours of held 
stars (J — Kg) and {H — Kg) versus Kg magnitudes along 
the reddening vector onto the reference RGB. The use of 
both colour-magnitude planes enabled us to obtain a mean 
extinction for each held and new constraints on the index 
of the near-infrared extinction power law a. We found that 
a value of a = 1.6 is inconsistent with the colours of inner 
Galactic stars, and by taking 47 Tuc as a reference for the 
RGB we hnd a = 1.9 ± 0.1. 

For Ag-band extinctions larger than 1.6 mag the 
2MASS (ATs, J — Kg) CMD provides extinction estimates 
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that are too low due to a selection effect from the J-band 
dropout of more distant sources. The 2MASS (it's, H — Ks) 
CMD suffers less from this bias. 

We reviewed near- and mid-infrared dereddened 
colour-colour relations of Mira stars and use them to test 
the quality of the extinction corrections for each SiO tar¬ 
get. Under the assumption that the SiO target stars are 
spatially distributed similar to the surrounding field stars, 
we corrected the photometric measurements of the SiO 
targets by adopting the median extinction of their sur¬ 
rounding field stars. The derredened colours of the SiO 
targets are not symmetrically distributed around the fidu¬ 
cial colour-colour line of known Mira stars, which suggests 
that for part of the SiO targets we may still be underes¬ 
timating the interstellar extinction by up to 15%. About 
50 SiO targets lie significantly in the “foreground” of the 
mean stellar distribution. 

Using the relation b etween mas s -loss rate and 
(As —15)o colour given by I.Teong et~ ] l)200,‘ljl . we esti¬ 
mated that most of the SiO targets have mass-loss rates 
in the range 10“^ to 10“^ Mq yr“^. 
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Appendix A: SAAO and 2MASS colours and 

magnitudes 

Transformation equations bet ween the colou rs and mag¬ 
nitudes measured in the SAAO ICarteil 1 19^ a nd 2MASS 
photo metric systems have been derived by ICarnent'^ 
(2^ using a list of mostly blue 9 4 photometric stan¬ 
dards. Figure 12 in ICarpentCT ll200lll shows that the dif¬ 
ferences between magnitudes and colours obtained with 
the two systems are smaller than 0.15 mag. 

Considering that Mira stars typically have a pulsa¬ 
tion amplitude in the near-infrared of 1-2 mag and that 
2MASS data are from a single-epoch observation ran¬ 
domly taken with respect to the stellar phase, the sys¬ 
tem transformations only have a secondary effect on the 
total colour and magnitude uncertainty, when compar¬ 
ing data from 2MASS with data taken with the SAAO 
telescope. However, since Mira stars are cold objects, we 
must exclude the idea that a combination of molecular 
bands and filter transmissions could generate a different 
colour transformation for this special class of objects. We 
therefore looked for 2MASS cou nterparts of the 104 oute r 
Bulge Mira stars monitored by I Whitelock et al.l lll994li . 
Mira stars are among the brightest objects detected in 
the Ks band such that identi fication o f their 2MASS coun¬ 
terparts is straightforward I Paper iJl (only three sources 
were excluded because they had no unique counterparts). 

Differences between t he mean magnitudes obtained 
with SAAO observations ilWhitelock et al.(ll99lh and the 
single-epoch 2MASS data have a dispersion of up to 0.8 
mag and the following mean differences: 


^s( 2 MASS) - ^(SAAO) = -0.15 ± 0.06 mag; 

(J - Ars)( 2 MASS) - {J - ^)(SAAO) = -0.14 ± 0.05 mag; 


{H - Ks)(^ 2 MASS) - {H - A')(saao) — -0.06 ± 0.03 mag. 

The colour-colour relations described in Sect. 0] using 
SAAO photometry also describe the 2MASS sequence of 
Mira stars well. 

Appendix B: IRAS and MSX filters 

Most of the past work was carried out using IRAS photom¬ 
etry; therefore, the currently available colour-colour rela¬ 
tions of Mira stars use mid-infrared data from the IRAS 
catalogue. A comparison of mid-infrared filters is therefore 
mandatory when translating old findings into new MSX 
and ISOGAL colours. 



C-D 

Fig. B.l. Colour-Colour diagrams of our SiO tar¬ 
gets. The IRAS 12^m magnitude is defined as [12] = 
—2.51ogFi2[Jy]/28.3. The continuous lines are our best 
fits. 


Figure Em shows the difference between MSX magni¬ 
tudes and IRAS 12 /rm magnitude for the SiO targets. 
Note that the D filter excludes the silicate feature around 
9.7 /im, while A and C filters include part of it (see Fig. 
E- Therefore the {D — [12]) colour shows a dependence on 
the (C — D) colour, which increases when the silicate fea¬ 
ture at 9.7 /rm starts to be self-absorbed. The {A — [12]) 
and (C — [12]) colours do not show any trend with the 
(C — D) colour. Due to both uncertainties of photometric 
measurements and source variability, scatter is large; how¬ 
ever we derived relations between the A, C, and D and 
the [12] magnitudes, as follow: 

A - [12] = 1.15 ±0.03 mag 
C - [12] = 0.05 ± 0.03 mag 
D - [12] = 0.08(±0.18) - 1.13(±0.05)(C' - D) mag. 




























